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Abstract We consider a (possibly) nonlinear interface 
problem in 2D and 3D, which is solved by use of vari- 
ous adaptive FEM-BEM coupling strategies, namely the 
Johnson-Nedelec coupling, the Biclak-MacCamy coup- 
ling, and Costabel's symmetric coupling. We provide a 
framework to prove that the continuous as well as the 
discrete Galcrkin solutions of these coupling methods ad- 
ditionally solve an appropriate operator equation with a 
Lipschitz continuous and strongly monotone operator. 
Therefore, the coupling formulations are well-defined, 
and the Galerkin solutions are quasi-optimal in the sense 
of a Cea-type lemma. For the respective Galerkin dis- 
cretizations with lowest-order polynomials, we provide 
reliable residual-based error estimators. Together with 
an estimator reduction property, we prove convergence of 
the adaptive FEM-BEM coupling methods. A key point 
for the proof of the estimator reduction are novel inverse- 
type estimates for the involved boundary integral oper- 
ators which are advertized. Numerical experiments con- 
clude the work and compare performance and cffcctivity 
of the three adaptive coupling procedures in the presence 
of generic singularities. 



1 Introduction 

1.1 Model problem 

Let Q C M. d (d = 2,3) be a bounded Lipschitz domain 
with polyhedral boundary T := df2 and normal vec- 
tor v. For given data (/, ito,<A)) € L 2 (Q) x H 1/2 (r) x 
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H 1 / 2 (J'), we consider the nonlinear interface problem 
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As usual, these equations are understood in the weak 
sense, i.e. we seek for a solution (u,u cxt ) G iJ x (i7) x 
H^ oc (n cxt ), with Hl c (Q cxt ) = {v\v e H X {K),K C 
12 ext compact}, and 12 ext := M d \77. It is well-known that 
problem (1) admits a unique solution in 3D. In 2D, the 
given data have to fulfill the compatibility condition 



(/.I) 



i,l)r = 0, 



(2) 



to ensure the right behaviour of the solution at infin- 
ity. Moreover, in 2D we assume diam(J?) < 1 to en- 
sure ellipticity of the simple-layer potential 23 defined 
below. The assumptions on the strongly monotone oper- 
ator 21 : R d -> R d will be given in Section 2.2. We denote 
by H^ 2 {r) the trace space of H 1 ^) and by H-^ 2 {r) 
its dual space. For simplicity, we will write u for the trace 
of a function u € if 1 (12), if the meaning is clear. 



1.2 Coupling of FEM and BEM 

Because of the presence of the unbounded exterior do- 
main l? ext in (lb), it is numerically attractive to rep- 
resent u cxt in terms of certain integral operators. This 
leads to a contribution on the coupling boundary r in- 
stead of the exterior domain i? oxt . For the interior do- 
main fl in (la), the (possible) nonlinearity of 21 as well 
as the (possible) inhomogeneity / ^ favours the use 
of a finite element approach. This led to the develop- 
ment of certain coupling procedures, and we focus on 
the Johnson-Nedclcc coupling [19], the Biclak-MacCamy 
coupling [8], and Costabel's symmetric coupling [13] in 
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the following. All of these approaches lead to a varia- 
tional formulation 

6(u,v) = £(v) for all v G H := H\f2) x H- 1/2 {r) (3) 

with unknown solution u£H, which is in a certain sense 
equivalent to (1). We equip % with the norm 

WAu = HW) + UW 2 H-vi{r) ( 4 ) 

for v = (v,tp) G H. Here, b{-, ■) is a continuous form on 
H x H which is linear in the second argument, and £(■) 
is a linear and continuous functional on T-L. The original 
works [8; 19; 13] focussed on linear 21 and hence bilinear 
&(•, •), and proved existence and uniqueness of the solu- 
tion u G T-L of (3). In the framework of the symmetric 
coupling well-poscdncss for nonlinear 21 has first been 
considered in the pioneering work [12]. 

For the Galerkin discretization, one considers a finite- 
dimensional and hence closed subspace He of H and seeks 
Ve E Hi such that 

b(Ui,V e ) =W) for all 6W<. (5) 

For linear 21, existence and uniqueness of the Galerkin 
solution U( G He for the symmetric coupling is already 
found in [13]. Moreover, Galerkin solutions are quasi- 
optimal in the sense of the Cea-type lemma 

||u - Ui\\n < C C 6a min ||u-V^||h, (6) 

where the constant Ccca > depends only on the geom- 
etry and on 21, but is independent of the given data, the 
continuous solution u, and the Galerkin solution U^. The 
analysis of the symmetric coupling has been generalized 
to nonlinear 21 in [12], but the proof required the under- 
lying mesh to be sufficiently fine, i.e. the maximal mesh- 
size had to be sufficiently small. Finally, for the nonsym- 
metric coupling strategies from [8; 19] and even linear 
problems, the analysis relied on the compactness of a 
certain integral operator & involved. However, this com- 
pactness restricted the coupling boundary to be smooth 
instead of piecewise polynomial. 

Only very recently, Sayas [23] proved that the Johnson- 
Nedelec coupling is equivalent to an elliptic problem, in- 
dependently of the compactness of .ft. For linear 21, more 
precisely the Yukawa or the Laplace equation, he thus 
derived that the variational formulation (3) as well as 
the discrete formulation (5) admit unique solutions and 
that the discrete solutions are quasi-optimal in the sense 
of (6). His analysis has been simplified by Steinbach [26]. 
For quite general linear 21, the latter work introduces a 
stabilized bilinear form 

b(u, v) := 6(u, v) + cr(u, v) (7) 

which is proved to be elliptic provided the smallest eigen- 
value of 21 is larger than 1/4. Up to some algebraic pre- 
/postprocessing, the solution u of (3) coincides with the 
solution u G H of 

b(u, v) = £(v) for all veH, (8) 



i.e. u = u + Uo. Steinbach thus proposed to approxi- 
mate the unique solution of (8) by some Galerkin solu- 
tion TJi G Hi and to obtain an approximation of u by 
XJe + Uo. One drawback of this method is, however, that 
the computation of the stabilization er(-, •) as well as of 
the (constant) offset Uo requires the (numerical) solution 
of an additional integral equation 5J</> C q = 1 • Firstly, this 
might lead to artificial error contributions from generic 
singularities of </> eq . Secondly, the first Strang lemma 
comes into play which imposes the assumption that the 
underlying (boundary) mesh is sufficiently fine. 

Finally, we mention the recent work [16], where for 
the (linear) Yukawa equation cllipticity of the bilinear 
form &(•, •) is proved for both the Johnson-Nedelec cou- 
pling as well as the Biclak-MacCamy coupling. 



1.3 A posteriori error estimation 

A posteriori error analysis aims to provide computable 
quantities gi which measure the Galerkin error ||u — 
Uf 11% from above (reliability) and below (efficiency). The 
local information provided by Qi can then be used to 
refine the mesh locally, where the Galerkin error ap- 
pears to be large. For the symmetric coupling, a pos- 
teriori error estimation was initiated by [12] for 2D and 
is well-established since then, cf. e.g. [10; 20; 27] and the 
references therein. To the best of our knowledge, only 
residual-based error estimators provide unconditional up- 
per bounds. On the other hand, for this type of estima- 
tors the lower bounds still require the mesh to be globally 
quasi-uniform although efficiency is also observed em- 
pirically on locally refined meshes [9]. For the Johnson- 
Nedelec coupling and the 2D Laplacian, different types of 
a posteriori error estimators have recently been provided 
and compared in [5]. 



1.4 Contributions of current work 

Adapting the results and proofs of [16; 23; 26], we present 
a framework which allows us to prove existence and unique- 
ness of the three coupling procedures for certain nonlin- 
ear 21. Roughly speaking, the idea is as follows: Each form 
b(u, v) on % which is linear in v, induces a nonlinear op- 
erator <B : % — > H*, where H* denotes the dual space of 
%. Then, the variational formulation (3) is rewritten in 
operator formulation 

<Bu = £. (9) 

For each coupling, we introduce an appropriate stabi- 
lization <r(-, ■) and consider the nonlinear operator *B in- 
duced by &(•,•) from (7). This is done in a way which 
ensures equivalence 

<Bu = £ ^ <Bu = £. (10) 
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Under appropriate assumptions on 21, the operator 03 
is Lipschitz continuous and strongly monotone (or: el- 
liptic). Therefore, the continuous operator formulation 
03 u = £ as well as its Galerkin formulation admit unique 
solutions u £ H rcsp. e He which also solve (3) 
resp. (5) and satisfy the Cea-type estimate (6). For the 
Johnson-Nedelec coupling and the Bielak-MacCamy cou- 
pling, our analysis requires that the cllipticity constant 
c n > of 21 is larger than 1/4, which reflects the same 
restriction as for the linear case in [26] . For the symmet- 
ric coupling, we avoid any restriction on c c n > 0. We thus 
obtain the same results as in [12], but without any re- 
striction on the mesh-size and with a much simpler proof. 
We stress that, unlike the approach of [26], the stabilized 
variant is only employed for theoretical reasons to guar- 
antee unique solvability of the non-stabilized equations. 

Finally, for lowest-order pieccwisc polynomials, wc 
derive residual-based a posteriori error estimators which 
provide reliable upper bounds for the respective Galerkin 
errors. For the Bielak-MacCamy coupling, we adapt the 
arguments from our recent preprint [4] to prove that the 
usual adaptive algorithm drives the residual error esti- 
mator to zero. 



1.5 Outline 

We start with a preliminary Section 2 which collects the 
precise assumptions on 21, the integral operators QJ, K, 
and 22J involved, as well as the notation used in the 
remainder of the paper. Section 3 then considers the 
Bielak-MacCamy coupling. We sketch the derivation of 
the coupling equations and prove existence and unique- 
ness of the continuous as well as of the Galerkin formu- 
lation as outlined above. Finally, we state and prove a 
residual-based a posteriori error estimator. In Section 4 
and Section 5, the same is done for the Johnson-Nedelec 
coupling as well as for Costabel's symmetric coupling. 
However, for the sake of brevity and since the proofs are 
very similar to that of the Bielak-MacCamy coupling, 
we only sketch the details. Emphasis is laid, however, 
on the fact that no restriction on the cllipticity constant 
c n > of 21 is imposed in the case of the symmetric cou- 
pling. Section 6 states the usual adaptive mesh-refining 
algorithm. Using the concept of estimator reduction and 
recent results of [4] , convergence of Ui to u is proved as 
£ — > oo, where I denotes the step counter of the adap- 
tive loop. A final Section 7 provides some numerical ex- 
periments. Emphasis is laid on the comparison of the 
three coupling procedures with respect to accuracy and 
computational time. Moreover, we numerically investi- 
gate the restriction c e ii > 1/4 in case of the Johnson- 
Nedelec and Bielak-MacCamy coupling. Finally, we see 
that the proposed adaptive schemes are much superior 
to the usual approach, where the mesh is only uniformly 
refined. 



2 Preliminaries 

2.1 Boundary integral operators 

Throughout, 8. denotes the double-layer potential with 
adjoint R> , 0J denotes the simple-layer potential, and 
2K the hypersingular operator. With the fundamental 
solution of the Laplacian 

_ J-^ log |z| forzeK 2 \{0}, 
[Z) -~ forzeK 3 \{0}, 

these integral operators formally read as follows, 

(Zl<p)(x) = J G(x-y)<p(y) dr y , 

{S*tp){x) = d v[y) G{x - y)(p(y) dr y , 

(Wip){x) = -d u(x) J d v ( y) G(x ~ y)f{y) dr y , 

for x E r and with d v ( y ) denoting the normal derivative 
at y G r. By continuous extension, we obtain bounded 
linear operators 

gjG£(^- 1/2 (r);ff 1/a (r)), 
&eL(H^(ry,H^(r)), 

^t e L(ff- 1 /2(r); J ff-V2 (r )) ) 

w e L(H 1/2 {ry,H~ 1/2 (r)). 

Finally, we stress the ellipticity of the simple-layer po- 
tential (<j>,%S<j))r > ||0]]^--i/2(n- Together with symme- 
try and continuity of 0J, this implies norm equivalence 
(</>,03</>)r — ||0||^_i/2(/y F° r further properties of the 
integral operators, the reader is referred to the litera- 
ture, e.g. the monographs [18; 21; 22; 25]. 

2.2 Strongly monotone operators 

An operator 03 : % — > %* is Lipschitz continuous pro- 
vided that there is a constant Cup > such that 

||®u-»v|| w . <C Lip ||u-v|| K (12) 

holds for all u, v £ H, where || • denotes the usual 
norm on the dual space H*. With (• , •} the duality brack- 
ets on H* x H, the operator 03 is strongly monotone 
provided that there is a constant C n > such that 

(03u-03v,u-v) >C u||u-v|& (13) 

holds for all u, v <E %. Wc refer to [29, Section 25.4] 
for the following standard results on strongly monotone 
operators: Under (12)— (13), the operator 03 is bijective, 
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and the inverse OS 1 is Lipschitz continuous with Lip- 
schitz constant 1/C e \\. Consequently, for every £ G H*, 
there is a unique u£H with 



(Q3u,v)=£(v) for all v e 



(14) 



and u depends Lipschitz continuously on £. Moreover, 
for every closed subspace He of H, there is a unique 
XJe G He such that 



(05 , V,) = £(Ve) for all V/ G He- 



(15) 



Finally, depends also Lipschitz continuously on £, 
and there holds the Cea-type quasi-optimality (6), where 

Ccoa = Chip/Cell. 

Remark 1 Provided that the discrete spaces He satisfy 



Remark 3 As far as existence and uniqueness of con- 
tinuous solution u and discrete solution XJe from (14)- 
(15) is concerned, our analysis only requires that the in- 
duced operator 03 is strictly monotone, i.e. 

(<8u-53v,u-v) > for all u,v G H with u ^ v, (20) 

i.e. (13) is replaced by (20). Then, the Browder-Minty 
theorem applies and, in particular, proves weak conver- 
gence \Je u in H as I — > oo under assumption (16). 
In this framework, however, the Cea-type estimate (6) 
cannot hold in general and a posteriori error estimates 
can hardly been derived. Therefore, we leave the details 
to the reader. However, we stress that (20) holds if the 
nonlinearity 21 satisfies 

o < (av« -avto, vv - vw) n (21) 

for all v, w G H (fi) with Vw ^ Vw instead of (19). 



He C He+i for all I > and {J^ He = H, 



(16) 



the quasi-optimality (6) implies convergence XJe — > u of 
the Galerkin solutions as t — > oo. In practice, the condi- 
tions (16) are satisfied if the underlying meshes are suc- 
cessively refined and the corresponding mesh-sizes tend 
to zero everywhere. 

To apply the framework of strongly monotone oper- 
ators to the FEM-BEM coupling formulations presented 
in Section 3, 4, and 5, we have to make some assump- 
tions on the coefficient function a : M. d — > M d . Firstly, we 
consider a to be Lipschitz continuous, i.e. there exists a 
constant cn p such that 



d . Integrating the square of (17), 



1% - M < ciipb - A (17) 

holds for all y, z G 
one obtains 

liavu - avHii a(n) < cf iP ll Vu - VHI!= W m 

for all v,w G i? 1 (J7). Secondly, we assume a to be 
strongly monotone in the following sense: There exists 
a constant c e ii > such that there holds 

\Vv - Vw\\ 2 L2 , n) < (aVw - 2lVro , Vv - Vw) n (19) 

for all v,w G H 1 ^). Here, denotes the L 2 (Q)- 

scalar product, i.e. (v,w)n = j Q vwdx. Similarly, we 
shall write {■ ,-)r for the L 2 (/ n )-scalar product which is 
extended by continuity to the duality bracket between 
H 1 ' 2 ^) and H" 1 ' 2 ^). 

Remark 2 We stress that conditions (18) and (19) are 
sufficient for solvability considerations, see Sections 3.2- 
3.3, 4-2~4-3, and 5.2-5.3. Anyhow, in our a posteriori 
analysis we need that a is pointwise Lipschitz continu- 
ous (17). 



Cell 



2.3 Discrete spaces 

In Sections 3-5, the model problem (1) is reformulated 
as variational equality (3) in the Hilbert space H := 
i? 1 (/2) xH~ 1 / 2 (r). For the respective discretizations, let 
7~e be a regular triangulation of ft and let Sf be a regular 
triangulation of T, where regularity is understood in the 
sense of Ciarlet. We approximate a function u G i? 1 (^2) 
by continuous, 7^-piecewise affine functions on Q. For a 
function <f> e H'^ir), we use £[ -piecewise constant 
functions, i.e. our discrete spaces read Hi := <S 1 (7i) x 
V°{£[) C H. 

Let £p denote the set of all interior faces, i.e. for 
E G £^ there exist unique elements T + ,T~ with E = 
T+ nr. We define the patch of E G £f by uj i>E := 
T + U T~ . Furthermore, we define the local mesh-width 
function he by 



he{r) := 



for r G 7*, 

fo r rebuff, 



where |-| denotes the volume resp. surface measure. A 
triangulation Te is called 7-shape regular, if there holds 

diam(T) d 



max ■ 

TeT e 



\T\ 



<7- 



Analogously we call £[ 7-shape regular, if 



max 

Eeel 



diam(E) 5 



E 



< 



(22) 



(23) 



for d = 3. For d = 2, the 7-shape regularity of £[ reads 



cr(£» ) := max 



I \E\ 



e' nE^d)} < 



7- 



(24) 



The definition of he and shape regularity implies equiv- 
alence diam(T) ~ he(r) ~ he(r') for all r,r' G 7^ resp. 
r, t' G £/" with t n t' ^ 0, where the hidden constants 
depend only on 7. 
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Remark 4 (i) We stress that Tit and £f are formally 
independent triangulations of fl and T, respectively. For 
the numerical implementation, however, we restrict to 
the case that £[ is the restriction Ti\r of Ti on the 
boundary, which indeed is a regular triangulation of T. 
In this case, we finally remark that "/-shape regularity of 
Te also implies j-shape regularity of £[ := Te\r ■ 
(ii) In 2D, the radiation condition (le) o/it GXt can also 
be adapted to u cxt (x) = a log |a;| + O(l) for x — > oo and 
fixed a G R. In this case, the compatibility condition (2) 
can be dropped. The analysis of the following sections 
still holds true for that case. 



3 Bielak-MacCamy coupling 

We can reformulate problem (1) with the help of the 
Bielak-MacCamy FEM-BEM coupling, which first ap- 
peared in [8]. This section is build up as follows: Firstly, 
we give a short sketch of the derivation of the Bielak- 
MacCamy coupling equations. Then, we investigate well- 
posedness of their continuous and discrete formulations. 
And last, we derive an residual-based error estimator for 
the Bielak-MacCamy coupling method. 



3.1 Derivation of Bielak-MacCamy coupling 

The first Green's formula for the interior part (la) reads 

(»Vu , Vv) a - (SIVu ■v,v)r = {f,v) n (25) 

for all v G -ff 1 (.!?). We plug in the jump condition (Id) 
for the normal derivative and obtain 

(%Vu,Vv)n-{Vu e * i -v,v) r = {f,v}n + (h,v) r . (26) 

For the exterior solution it oxt of (lb), we make an indirect 
potential ansatz with the simple-layer potential 



u ext = 030 in n e 



where the integral operator 03 is defined as 03, but is now 
evaluated in J7 GXt instead of T. We stress that the density 
(j> G H~^^{T) is unknown. Then, we use properties of 
the simple-layer potential operator: Firstly, by use of the 
continuity of the simple-layer potential in R d , i.e. 030 = 
030 on T, and the trace jump condition (lc), we see 



-u + 030 = — uq on T. 



(28) 



Secondly, we use the jump condition of the exterior conor- 
mal derivative of the simple-layer potential to see 



Vu GXt • v = 9 GXt O30 = -(A - £t>. 



Plugging the last equation into (26) and supplementing 
the system with the variational formulation of (28), we 



end up with the variational formulation of the Bielak- 
MacCamy coupling: Find u = (u, 0) G H such that 



{%SJu,Vv) n + {{\-&)4>,v)r 
-<V>,u)r + (^,Sty)r 



{f,v)n + (<f>o,v} r , 
-(iP,u )r, (30) 



holds for all v = (v, G H. 

From now on, let Xe be a closed subspace of H 1 (H) 
and yi be a closed subspace of 77 -1 / 2 (-T). We define 
He := Xi x ye. Note that the entire space H = He is a 
valid choice, and hence the following analysis applies to 
both, the continuous formulation (30) and the Galerkin 
discretization. In the latter case, u G H in (30) is re- 
placed by G He, and v G H is replaced by arbitrary 
V, G H e . 



3.2 Stabilization 

We define the linear form 6bmc 
u = (u, 0), v = (v, if)) G H by 



HxH 



for any 



6 bmc (u, v) := (!*Vu , Vv) n + <(| - #)0 , v) r 



(31) 



Note that 6bmc( - , •) is only linear in the second argument. 
Furthermore, we define linear functionals a\ and a-i on 
H x {n) and H" 1 I' 1 {F) by 



a2W0 := (V , .-'"o)r 



(32) 
(33) 



for all v = (v, iji) G H. For 02, we also use the notation 
a2(V0 = (V 1 j a 2) r- With these definitions, the continuous 
formulation of the Bielak-MacCamy coupling is equiva- 
lently written as follows: Find u G H such that 



&bmc(u, v) = £(v) := oi(u) + a 2 (jp) 



(34) 



for all v G H. Moreover, the Galerkin formulation of 
problem (30) reads: Find XJe G He such that 



(27) b hmc (V i ,V e ) = Z(V i 



(35) 



holds for all = {V t , &i) G He- 

Throughout the remainder of this section, we need 
the following assumption. 

Assumption 5 There is a fixed function £ G H^gN ^ 
with (£,l) r ^0. 

Remark 6 The discrete space ye = V°(£f), introduced 
in Section 2.3, fulfills Assumption 5 with £ = 1. 

Now, we try to show ellipticity of a linear form which 
is equivalent to 6bmc(', •)• Firstly, note that we have to 
take care of the fact that &bmc(', •) is not elliptic since 



(29) bbmc((l, 0), (1, 0)) = (avi , Vl)n = 0. 



(36) 



Therefore, we introduce a new linear form &bmc(', •) which 
is equivalent to 6bmc( - 7 •)• 
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Theorem 7 With £ from Assumption 5, the linear form 



frbmc(U£, V e ) := 6 b mc(Uf, Vi 



(37) 



is equivalent to the linear form &bmc( - j •) w& f&e following 
sense: The pair XJe = (Ue,<Pe) G solves problem (35) 
i/ and on??/ i/ li solves 

&bmc(U^, V<) - £(V/) + (C , a 2 ) r (£ , QJtf', - V e ) r (38) 
/or V^^^eHf. 

Proof Step 1. Let U £ = (Ut,$t) be a solution of (34). 
Testing with (Vfc.lfc) = (0,0 G we see b b mc(U^ (0,0) 
= 02(C)- With the definition of 6bmc(', Oi we infer 

= 6 bmc (U,, (0,0) - 02 (0 = (t,mt-U t -02)r- 

Hence, (£,2J<2> f - U e - a 2 )r{€ ,W&e - V e ) r = for all 
~Ve & He- Clearly, this is equivalent to 

(£ , m t - U e )r{i , m t - Vf) r = (C , o 2 )r(C , «W< - Vt)r 

for all V £ = (Vt,&t) G Therefore, = (U t ,$ t ) G 
"H^ also solves problem (38). 

Step 2. For the converse implication, let XJe = (Ue, $g) G 
He solve (38). The choice of V e = (0,0 in (38) gives 

a2(O + (Oa2)r(O s »Or = &bmc((C/f,^),(0,O) 

= (£ , m t - U t )r + « , - ^)r(C , WQr, 

which is equivalent to 

(£ , - U t - a 2 ) r (l + (£ , !D0r) = 0. 

Since 03 is ff _1 ' 2 (P)-elliptic, the last equation implies 
(f , <U<^ -U t - a 2 ) = 0. Wc thus infer 

u bmc (Ve,Ve) 

= (f,9J^-D/)r^,aW>i-Vi>r 
= (Oa2)r(0^-^)r 

which, together with (37) and (38), proves that Ve = 
(Ue,<Pe) is also a solution of problem (34). □ 



3.3 Existence and uniqueness of solutions 

The linear forms Ob mc (-, ■) and 6 bmc (-, •) induce operators 
®bmc, Sbmc :H^H* by 



(SbmcU , v) := 6 b mc(u, v) for all u, v G H, 
(QSbmcU , v) := 6bmc(u, v) for all u, v G H. 
The main result of this section reads as follows: 



(39) 



Theorem 8 Under Assumption 5 and provided that the 
ellipticity constant c n of 21 fulfills c c \\ > 1 /4, the opera- 
tor Q5bmc is strongly monotone and Lipschitz continuous. 



The proof requires the following lemma which is proved 
by means of a Rcllich compactness argument. 

Lemma 9 For u = (u, <fi) G H, let 



||u|| 2 := ||V U » 2 



L2(n) + (j>Mr + |(e,9ty-u) r | 2 , (40) 



where £ is provided by Assumption 5. Then, ||| • ||| defines 
an equivalent norm on H. 

Proof Clearly, there holds |||u||| < ||u||^ for all u = (u, cf>) 
G H. To see the converse estimate, we argue by con- 
tradiction and assume that ||u n ||^ > «.|||u„||| for certain 
u„ = (tin, <f> n ) and all n G N. We define v„ = (v n , ip n ) by 



\\Un\\n 

and obtain \\(v n , tp n )\\n = 1 as well as |||(i;„, ip n )\\ < 1/n. 
By definition of ||| • ||| and ellipticity of 53, this implies 
ip n — >• G H-^ 2 (r) and Vu„ -> G i 2 (!2). More- 
over, by extracting a subsequence, we may assume that 
(v n ,i/} n ) (i>, VO in H. Clearly, this implies ip = and 
w„ — 1 v G H 1 (il), whence i>„ — > v G L 2 (Q). Moreover, 
weak lower semi-continuity of ||| ■ ||| implies |||(t>,'!/>)||| = 0; 
whence Vw = and |(£,v)r| = 0. From the choice 
of £ and since v is constant, we infer v = and thus 
»„^0e H 1 ^). Altogether, v„ = (v n ,tp n ) -> G % 
contradicts ||v n ||^ = 1. □ 

The following proof of Theorem 8 is very much influ- 
enced by the investigations of [23; 26]. We recall some ba- 
sic facts on the boundary integral operators, cf. e.g. [25, 

Chapter 6]: For given x £ #~ 1/2 (-0, let u * = ^X- Then 
there holds 

• d v u* = (\ + 8J)x, 

• (d v u* ,v) r — (V«, ,Vv)n for all v G H X (Q), 

• (X , «x)r = ([cU*] , 53 X )r = || Vu. ||| 2(Rli) , 

where the last equation is only valid for d = 2, if x S 
H~ 1/2 (r) = {ipe H-V 2 (r) I (V>, l)r = 0}. For arbitrary 
X G i? _1//2 (/ n ), we therefore introduce the splitting 

X = X*+Xeq, (41) 

with x* G F~ 1/2 (r) and x cq = (x,l)r </> Gq .Here, oq = 
53 _1 1/(2J _1 1 , l)r is the so-called equilibrium density 

— 1/2 

(or: natural density). Note that x* £ ^* (^)- More- 
over, with the commutativity relation ^qj-i = QJ -1 ^ 
and the equality (h + R)l = 0, we infer 

3 (42) 

for all u G H 1/2 (T). Together with the splitting (41), 
this proves 

((l + tf) X ,v) r = ((± + tf)x*,v) r . (43) 
Finally, there holds 

(X* , 2?Xc q )r = (X* , l)r(x , l)/^/^ 1 ! , l)r = 
and therefore 

( X , X X )r = (X* , %X*)r + {Xo q , »Xc q )r. (44) 



((i + S^- 1 ! , v) r = (23- X (| + ^)1 , «)j 
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Proof (of Theorem 8) Lipschitz continuity of OSbmc sim- 
ply follows from the Lipschitz continuity of 21 and the 
continuity of the boundary integral operators. 

It thus only remains to show ellipticity of 2$bmc- Let 
u= (u,(f>),v= (v,4>) e H. Then 

(SbmcU - S bmc V , U - v) 

+ ((i -$)(<j>-ip),u- v) r -{<f>-i/>,u-v) r 
+ (0 - V , ST(0 - V0)r + |<f , !0(^ - ^) - (« - «))r| 2 
=:7i+/ 2 + / 3 +/4 + l5. (45) 

Below, we show 

/1+/2+/3+/4 > ||V«-Vt;||i a(n) + (0-V,2J(0-V))r. 
With Lemma 9 and the definition of ^5, this implies 

(QSbmcU - QSbmcV , U - v) > |||u - v||| 2 

and thus concludes the proof. 

Step 1. To abbreviate the notation, we write w = (w, x) 
= u — v. The term Ii is estimated by strong monotonic- 
ity (19) of 21, 

(2lVu-2lV«,Vu;)fi > c e ii|| Vw\\ 2 L2{n) . (46) 

Step 2. With the splitting (41) of \ and u* = 23x*i tnc 
terms J 2 + ^3 can be estimated by 

h + h = -{{\+ f#)x , ™)r = - ((3 + £ f )x* , ^>r 

> -||Vu*|| L 2 (r2) ||Vw|| L 2 (fi) 

> -||V«*||i2( K *)||Vl0||£2(fl). 

Step 3. Wc recall Young's inequality: For arbitrary a, b G 
E, 5 > there holds ab < § a 2 + ^-b 2 . We infer 

- ||Vu*||z,2(R«i)||Vlu||i2( fi ) 

> -|l|V^|||, (r2} - A_i|| V ^||i 2(Rd) . 

We combine the second term with I4 and see 

-^||Vu*|| 2 2(Rd) + (x,2Jx)r 

= -^<X* , %X*)r + (X* , W X *)r + (x cq , WXc q )r 
= (1 - i 3 i )(x* ,%X*)r + (Xc q ,5JXc q )r 

> (l-^2 i )(x,»x)r, (48) 
Step 4. We combine (45)-(48) and obtain 

(SbmcU - §bmcV , w) > (c e ll - f) || VtO 1 1 ^ 2 {n) 

+ (l- i 3 i )(x,5Jx)r 
+ |(e,2Jx-w;)r| 2 . 



We have assumed that c c n > 1/4. Hence, there exists 
some 8 > with 1/4 < 8/2 < c c \\. Furthermore, such a 

(5 implies c c ii — § > as well as 1 — > 0. We define 

C c n := min{c c n — |, 1 — } an d end up with 

(<B b mcU - SbmcV,w) > C eU |||w||| 2 . 

With Lemma 9, this proves ellipticity of *Bbmc- □ 
Remark 10 (i) In the case d = 3, there holds 
\\Vu4l 2m = ( X ,%x)r 

for all x € i/ _1 / 2 (-T) anc? 11, = 2Jx- Then, the proof 
of Theorem 8 simplifies, because the splitting (41) is not 
needed and one may simply choose x — X* ■ 
(ii) The assumption c c n > 1/4 is sufficient, but may not 
be necessary. Numerical experiments for a linear opera- 
tor 21 have shown that the bound 1/4 is not sharp, i.e. 
solvability seems to be given also for < c c n < 1/4. 

Finally, we may apply the standard results from the 
theory on strongly monotone operators, see Section 2.2, 
to prove in conjunction with Theorem 7 the following 
corollary. 

Corollary 11 Under the assumptions of Theorem 8, the 
Bielak-MacCamy coupling (34) admits a unique solution 
u £ H. Moreover, Galerkin approximations TJg G Hi 
of (35) are quasi-optimal in the sense of (6). 

Proof We define <8 := !Bb mc , and let £ be the right- 
hand side of (38). According to the main theorem on 
strongly monotone operators, the operator equation (14) 
and its Galerkin discretization (15) admit unique so- 
lutions u G H and XJe G He- Moreover, these satisfy 
the quasi-optimality (6). Finally, Theorem 7 proves that 
ueWis the unique solution of (34) , and \Jg G He is the 
unique solution of (35). □ 



3.4 Residual-based error estimator 

Our aim is to derive a reliable residual-based error es- 
timator for the Bielak-MacCamy coupling in the same 
manner as in e.g [12] or [5]. 

Let [2lVlfy -vWe denote the jump of 2lV£/g ■ v over the 
interior face E £ £p . We assume additional regularity 
4>o £ L> 2 {r) and uq G H 1 ^) from now on. 

Theorem 12 Suppose that u G H is the unique solution 
of the Bielak-MacCamy coupling (30), and \Jg G He := 
«5 1 (7]?) x V°(£f) is its Galerkin approximation. Then, 
there holds 

CTefllu - U.H 2 , < p 2 := MT) 2 + 

TeTe E&SpUSf 



The volume contributions read 

P z(T) 2 = h 2 T \\f\\h iT ) forTeTe, (49) 
p e (E) 2 =h E \\[WUfv}\\ 2 L2{E) }orEe£l\ (50) 

whereas the boundary contributions read 
Pi (E) 2 = h E \\<j> + (Jit - - %VU e ■ v\\ 2 L2{E) 
+ h E \\Vr(U e -u ~ < 8$t)\\ 2 L 2 (E) 

for E G £[ ■ The constant C Te \ > depends only on Q, 
r and the "{-shape regularity of Tit ^ n d £[ ■ The sym- 
bol Vr(') denotes the surface gradient (resp. arclength 
derivative for d = 2). 

Proof Recall the definitions (32)-(34) and (39) of 93 bmc , 

£, and %$bmc- Problem (30) for the exact solution and 
its Galerkin approximation arc cquivalently written as 

(Q3 bmc u,v) = £(v) forallvGH, 

and 

(*B bmc U, , V,) = £(V<) for all V, G U t . 

Wc stress that (<8 b mcU - QSbmcU^ , u - U £ ) - (*8bm C u - 
Q3 bmc U, , u - U,) = \(^V3(<j)-$ e )-(u-Ue))r\ 2 = 0, 
which follows from the second equation of the Bielak- 
MacCamy coupling, cf. (30). With ellipticity of 93 bmc , 
we obtain 

[u - U/|& < (S bmc u - 5 bmc U £ , u - U £ ) 

= (QSbmcU - S B b mcU £ , U - U e ) 

< HQSbmcU - S bm cU,||„, || U - V e \\ n - 
This leads us to 

||u - U e \\ n < ||Q5bmcU - B^U/I^. 

|£(v)-<93 bni cU^v}| 



sup 

vew\{o} 

sup 

v£«\{0} 



|£(v-V/)-(»bmcU <I V-V 4 )| 



l v llw 



for arbitrary G We choose Vi = 0) G 
where Zi ■ H 1 (f2) — >• 5 1 (7i) denotes a Clement-type 
quasi-interpolation operator, which satisfies a local first- 
order approximation property 

||w - &HU 2 (t) $ diam(T)||Vto||£2( WT ) (51) 
and local iJ 1 -stability 

||V(u;-^iy)|U=(r) < I|V«;|| L 2 (wt) (52) 

for all io G H 1 ^) and T G Tt- Here, w T = U{ r ' G 
7e | T' n T ^ 0} denotes the patch of an element T eT e . 
An example for such an operator Zt is the Clement op- 
erator [1] or the Scott-Zhang projection [24]. Note that 



the constants in the estimates (51)-(52) only depend on 
fl and the 7-shape regularity of Ti- An immediate con- 
sequence of these properties and the trace inequality is 
that 



\\w - Ztw\\ L 2 (E) < &ia,rt\(E) 1/2 \\w\\ H i {u]T) , 



(53) 



where E is a face of the element T G Ti with E C T. 
Note that 



£(v- V,)- (<B bmc U,,v- V,) 
= {f,v-Ztv)n - (aVtf/,V(t> 
+ (0o - {\ - , v - 2tv)r 

+ (ip,Uz-u -Z}$t} r . 



(54) 



The first term on the right-hand side is estimated by use 
of Cauchy inequalities and (51). This gives 

(f,v- 3iv)n = E ^ ' w ~ Ztv) T 
TeTt 

< E ll/IU 2 (r)l|w-3HU 2 (T) 
TeT ( 



< E diam(T)||/|| i2(T) ||V«|[ L 2 



TeTi 



<(Eiimi! 2( t)) 1/2 (Eii v -ii^)) 

TGTe TeTe 
^ ||^/IU 2 (/2)||V«||i2 (fi ). 

Piecewise integration by parts of the second term on the 
right-hand side of (54) yields 



1/2 



TeT e 

= E (~ div 2lVC/f , u - &u)t + (2tVf/^ ■ 1/ , u — & v) dT 



TeTe 



= (%VUfV,v-3iv)r+ E ([^ u M,v-3tv) E , 



Ee£ 



since div 21VL^ vanishes clcmentwise because of 21VL^ G 
V°(Te). The second and third term on the right-hand side 
of (54) can thus be estimated by 

(0o - (§ - £t)^ - aVtf/ ■ f , « - 

iit)^ - 2tW £ ■ H|z,2 (B) ||v - &«]|z, 2 (e:) 

- E W\^ vu f ^Wl 2 (e)\\v -3M\li(e) 
Eeep 

--: Ji + J-2 



EGSP 

< E 11^0 -(5 

Beef 
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For each boundary face E we fix the unique element Te 
with E C Te and infer by use of (53) 



Ji < diam(.E) 



1/2 



< ||/i] /2 (0o - (I - £ f )^ - 3lVtf< • v)\\ L 2 ( r)\\v\\m(n). 

For an interior face E we fix some element Te with _B C 
Tg and estimate the term J 2 analogously by 



J 2 < 



E 



\hl /2 [WUt-v]\\hwy%\\mw 



To estimate the fourth term in (54) we observe 
(l,U e -u -Z3$ e ) E = 

for any E G 1?/", which follows from the second equa- 
tion of (30) tested with the characteristic function of E, 
which belongs to P°(£[). Now, [11, Corollary 4.2] can 
be applied and proves 

||f^-it -QJ^||^ 1/2(r) 

<Cio C J] cnam(£)||Vr(^-u -23^)||£ 2(£) , 
Ee£[ 

where C\ oc > depends only on r and the 7-shape reg- 
ularity of £[ . This leads us to 

(V; ,U t -u - %$e)r 

< \\lp\\ H -l/2 {r) \\Ul - U - W$l\\ H l/2 {r) 

< \\1>\\H-vnr)\\h] /2 Vr(Pi -u - m e )\\ LH r)- 

In 2D, the same estimate can also be obtained by use 
of the continuity of U t - u - 53<Z> £ G H 1 ^), cf. [12]. 
Altogether, we have 

|£( V - V t ) - (SbmcUf , v - V e }\ /\\y\\ H 

< l|WIU^) + ||^ 1/2 [2tV^-z,]|| i2(U£f) 

+ ||^ 1/2 (0 O - (I - - SIVDi • ^||^(r) 

+ ||^ 1/2 Vr(C/,- U o-53^)|| L 2 (r) , 

which concludes the proof. □ 



4 Johnson-Nedelec coupling 

In this section, we present the Johnson-Nedelec coupling, 
which first appeared in [19]. As in the previous section, 
we state the continuous and discrete formulation of this 
method and discuss existence and uniqueness of the cor- 
responding solutions. Finally, we provide a reliable resi- 
dual-based error estimator. 



4.1 Derivation of Johnson-Nedelec coupling 

Unlike the Bielak-MacCamy coupling, we represent the 
exterior solution u cxt by use of the third Green's identity 
in the exterior domain J? cxt , 



£u ext - 230, 



(55) 



where we define = <9„u cxt . As above 23 and M. are de- 
fined as 23 and but are now evaluated in J? ext instead 
of r. Taking the trace in (55), we see 

u e* = (j| + !) lt «a_gty. (56) 

Using the trace jump condition (lc), we obtain 

(i -£)u + 230 = (| -&)u . (57) 

Together with (26), the variational formulation of the 
latter equation provides the Johnson-Nedelec coupling: 
Find u = (u, 0) G H such that 



(58) 



(2lVu, Vv) u - (<f>,v) r = (f,v) a + (0 o ,v)i 
holds for all v = (v, -0) G H. 
4.2 Stabilization 



For the Johnson-Nedelec equations, we can apply similar 
techniques and derive similar results as for the Bielak- 
MacCamy coupling, see Section 3. For the sake of com- 
pleteness, we state these results in the following. 
We define the linear form 6j n (-, ■) for u, v G H by 

6jn(u,v) := (mu,Vv} n - {<f>,v) r 

+ (l& > (f-J*)u + aty) r . 

Furthermore, we define linear Junctionals a\ and a 2 by 

a i( v ) '■= (f,v)n + (<p ,v)r, 
a 2 (V>) := (V>) (§ - &)v>o)r 

for all (v, xp) G H. Then, problem (58) can be reformu- 
lated: Find u G H such that 



bj n (u, v) = £(v) := a±(v) + a 2 (ip) 



(59) 



holds for all v = (v, ip) G H. Moreover, the Galcrkin 
discretization of (59) reads: Find G He such that 



b jn (XJ t ,V t ) =£(V, 



(60) 



holds for all G He, where He = Xg X ye is a closed 
subspace of H. Similarly to Theorem 7, one proves the 
following result: 
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Theorem 13 With £ of Assumption 5, the linear form 4.4 Residual-based error estimator 



6 jn (u,v) := 6j n (u,v) (61) 
+ <£,(!- R)u + 2J0)r{£ , (| - £)« + 2J<A)r, 

is equivalent to the linear form &j n (- , •) m i/ie following 
sense: The pair TJg = (Ue,<I>e) € 7^ solves problem (60) 
i/ and cmiy i/ zi solves 



(62) 



MU*V*) = oi(^) + a3(#<) 

+ (e,« 2 )r(£,(^^)^ + QJ^)r 

for allV i = {V e ,&t) eH*. □ 

4.3 Existence and uniqueness of solutions 

We stress that there is a close link between the Johnson- 
Nedelec and the Bielak-MacCamy coupling, since 

6j n (u, u) = 6 b mc(u, u) for all ueM. 

This indicates that the analytical techniques to prove 
cllipticity of the two coupling methods are similar. 

The stabilized bilinear form 6j n (-, •) of (61) induces a 

nonlinear operator 23j n : % — > %* by 



&jn(u,v) 



(63) 



for all u,v G H. The following theorem states strong 
monotonicity of Q3j n under the same assumptions as for 
Theorem 8. Instead of Lemma 9, we need the following 
result: Under Assumption 5, the definition 



l Vu ll|2(fi) 



+ \(t l a-&)u+vs<t>) 1 



(64) 



for u = (u, 0) G 'H provides an equivalent norm on W. 
The proof is achieved by a Rellich compactness argument 
as in the proof of Lemma 9. 

Theorem 14 Under Assumption 5 and provided that 
the ellipticity constant c e \\ of 21 fulfills c c \\ > 1 /4, the 

operator Q3j n is strongly monotone and Lipschitz contin- 
uous. □ 

As above, Theorem 13 and standard theory on strong- 
ly monotone operators now prove the following corollary. 

Corollary 15 Under the assumptions of Theorem 14, 
the Johnson-Nedelec coupling (59) and its Galerkin dis- 
cretization (60) admit unique solutions u£H resp. G 
%£. Moreover, there holds quasi- optimality (6). □ 



We recall a residual-based error estimator for the Johnson- 
Nedelec coupling. The proof of the following result can 
be achieved with similar techniques as in the proof of 
Theorem 12. The linear case for d = 2 is found in [5]. 

Theorem 16 Suppose that u G T~L is the unique solution 
of the Johnson-Nedelec coupling (58) and G %i = 
S 1 (Te)xV°(£[) is its Galerkin approximation (60). Then, 
there holds reliability in the sense of 



C7el 2 ||u-U4«<% 2 



m(Tf 



E 

Eeepue' 



The volume contributions for % and £^ are the same as 
above, cf. (49) and (50), but are denoted by rji{T) resp. 
r)t(E). The boundary contributions read 

m {Ef = h E \\<h +*/ - 21V(7, ■ v\\l* {E) 

+ M|Vr((| - £)(uo - U t ) - W$i)\\h {E) 

for E G S[ . The constant C rc i > depends only on Q, 
F and the 7 -shape regularity of Tit an d £f ■ d 



5 Costabel's symmetric coupling 

In this section, we treat the symmetric FEM-BEM cou- 
pling method, which first appeared in [13]. First, we 
consider the continuous and discrete formulation. Af- 
terwards, we investigate well-posedness of the coupling 
equations, where we give a simpler proof of unique solv- 
ability than in the pioneering work [12]. Finally, we recall 
a residual-based error estimator from the latter work. 



5.1 Derivation of symmetric coupling 

We start from the Johnson-Nedelec coupling (58) and 
modify the first equation: With the ansatz (55) for u ext , 
we use the second Calderon identity 

= ^ x V xt = -(^ t -i)^-2Hu ext . (65) 

The trace jump condition (lc) eliminates it oxt and gives 



= -(JF - i)0-2Di 



2Bu n 



(66) 



We plug this identity into the first equation of (58) and 
move %Buq to the right-hand side. Altogether, the vari- 
ational formulation of the symmetric coupling reads as 
follows: Find u = (u, <f>) G % such that 



(2tVu , Vv) n + ((& ~±)<j),v)r + (Wu,v}r 

= (f,v)n + (<po + Wu ,v}r, (67) 

*)«)r + (i> 1 %<t>)r = W , (| - $)u )r 

hold for all v = (v, ip) G H. Note that the second equa- 
tion of (67) is just the same as for the Johnson-Nedelec 
coupling (58). 
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5.2 Stabilization 



5.3 Existence and uniqueness of solutions 



Unique solvability for the symmetric coupling (67) as 
well as for its discretization can be found in [12] for 
our nonlinear model problem (1). Nevertheless, here we 
present a much simplified proof of this result. We pro- 
ceed as before and define the linear form b sym (•, •) for all 
u,veWby 

6 sym (u,v) := (%.Vu,Vv) n + ({& -%)(/>, v) r 

+ (Wu,v) r + {ip,{j-Si)u + Z]4>)r. 



Moreover, we define linear functionals a\ and a 2 by 

ai(v) := (f ,v)n + (4>o + Wu a ,v) r , 
02 WO := (V>, (5 - &)u )r- 

Then, problem (67) can equivalently be stated as: Find 
u£W such that 



6 sym (u, v) = £(v) := oi(u) + a 2 (VO 



(69) 



holds for all v = (v, ijj) & T~L- For the discretization, let 
~Hi = X(_ x be a closed subspace of %. The Galcrkin 
discretization of (69) then reads: Find XJe G He such 
that 



b Bym {V t ,V e ) =£(V e 



(70) 



holds for all G He- 

The following theorem states equivalence of the lin- 
ear form 6 sym (-, •) to some new stabilized form 6 8ym ( - , •) 
which will turn out to be elliptic. The proof is achieved 
by similar techniques as used for proving Theorem 7 and 
is thus omitted. 

Theorem 17 With £ of Assumption 5, the linear form 



i( U , v ) 



i( u > v ) 



(71) 



The stabilized bilinear form 6 S ym(', •) from (71) induces 
a nonlinear operator Q3 sym : H — > H* by 



(93 sym u,v) := 6 sym (u,v) 



(73) 



for all u, v G H. The following theorem states strong 
monotonicity of Q5 sy ni without any further restriction on 
the ellipticity constant c e n > of the nonlinearity 21. 

Theorem 18 Under Assumption 5, the operator Q5 S y m 
is strongly monotone and Lipschitz continuous. 

Proof Lipschitz continuity of Q5 sym simply follows from 
the Lipschitz continuity of 21 and the continuity of the 
boundary integral operators. 

It thus only remains to show ellipticity of 23 sym . We 
write w = u — v = (w, x). Recall that the hypcrsingular 
operator 2U is positive semi-definite. With the ellipticity 
constant c e n > of 21, we see 



(03 



(74) 



2 



is equivalent to the linear form 6 sym (-, •) in the following 
sense: The pair TJe = {Ue^e) € Tie solves problem (70) 
if and only if it solves 



sym" -^sym v ; vv / 

= (%Vu-X.Vv,Vw)n + ((S# - \)xMr 
+ {Ww , w) r -<xM~ \) w )r + (x , %x)r 
+ \(t,V3 X +(k-&)w)r\ 2 

> c cll ||V^||| 2(r2) + <X,^X>J 

+ |(e,5Jx+(|-^V)r| 
>min{ CcI1 ,l}(||V W ||2 2(fi) + ( x ,QJ x ) r 

+ |(C,2Jx+(i-^Mr| 2 ). 

Norm equivalence for the norm ||| • ||| defined in (64) thus 

concludes the proof of ellipticity of the operator Q3 sym . 

□ 

As above, Theorem 17 and standard theory on strong- 
ly monotone operators prove the following corollary. 

Corollary 19 Under Assumption 5, the symmetric cou- 
pling (69) and its Galerkin discretization (70) admit 
unique solutions u G % andXJe G He, respectively. More- 
over, there holds the quasi- optimality (6). □ 

Remark 20 (i) Note that there is no restriction on the 
monotonicity constant c c \\ > of 21 as is the case for the 
Bielak-MacCamy and the Johnson-Nedelec coupling, 
(ii) In contrast to [12], our analysis does not need the 
mesh-size he to be sufficiently small to prove ellipticity 
of the coupling equations. 



b sym (Ue,Ve) = £(V e ) 
for all V e = (V e ,&t) G He- 



(72) 



□ 



5.4 Residual-based error estimator 

The proof of the following result can be found in [12] for 
d = 2 and is achieved by similar techniques as in the 
proof of Theorem 12. It is therefore omitted. 
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Theorem 21 Suppose that u € Ti is the unique solution 
of the symmetric coupling (69) and Ui G "He — 5 1 (J7) x 
V°(£[) is its Galerkin approximation (70). Then, there 
holds reliability in the sense of 



C, 



rel ll U - U dl« < n 



TGTt 



E 

Eeepus' 



Assumption 23 (Mesh- refinement) Let (%,£[) be 
a sequence of regular triangulations, where Ti+i is ob- 
tained from l~i and £f +1 from £[ by local mesh-refinement. 
Then, there are (-independent constants < 7 < 00 and 
< q < 1 such that the following holds 



The volume contributions for Ti and £f are the same as 
above, cf. (49) and (50), but are denoted by ni(T) resp. 
[ii(E). The boundary contributions read 

ME) 2 =h E \\<t>o - aW* • v + W(u - U t ) 



The triangulations (T},£j ) are j-shape regular. 

't 



A refined element T G T~i resp. face E G £[ is the 



- (£t - ^)$e\\l HE) £[ satisfy \E'\ < q\E\ 



union of its sons T' G Te+i resp. E' G £[ +1 - 
The sons T' G Te.+\ of a refined element T G Te sat- 
isfy \T'\ < q\T\. The sons E' G £p +1 of a refined 
interior face E G £f satisfy \E'\ < q\E\. 
The sons E' G £f +1 of a refined boundary face E G 

□ 



+h E \\Vr{(h ~ £)(«o - U t ) - 



li 2 (B) 



for E G Sf . The constant C rc \ > depends only on Q, 
T, and the j-shape regularity of 7} and £[ . □ 



6 Convergence of adaptive scheme 

In this section, we state an adaptive algorithm for the 
three different coupling methods and prove convergence 
of the adaptive scheme for the Bielak-MacCamy cou- 
pling. 



6.1 Adaptive algorithm & Mesh-refinement 

In the following, we fix one particular coupling and let 
Ue G Tie := S x {Ti) x V°{£[) denote the correspond- 
ing Galerkin solution, i.e. solves either problem (35), 
(60), or (70). By Q, we denote the corresponding residual- 
based error estimator, cf. Sections 3.4, 4.4, and 5.4. Un- 
der these assumptions, the standard adaptive scheme 
reads as follows: 

Algorithm 22 Input: Initial triangulation (To,£F), 
counter £ := 0, and adaptivity parameter < 8 < 1. 

(i) Compute discrete solution G He- 
(ii) Compute local refinement indicators (i(t) for all r G 
TeU£pU£[. 

(Hi) Find (minimal) set Aie Q Te Uf^ U£f such that the 
error estimator Q fulfills the Dorfler marking 



E 



UE) 2 . (75) 



TeM e nT t 



EEMen(EfUSf) 



(iv) Obtain new mesh (7~e + i,£f +1 ) by refining at least all 

elements of the set M.e- 
(v) Increase counter £ by one and goto (i). 

Output: Sequence of Galerkin solutions (U^)^ g N o-nd 
sequence of error estimators (Q)eeN - □ 

To prove convergence of the adaptive algorithm, we 
need certain assumptions on the mesh-refining strategy 
used in step (iv). 



We note that Assumption 23 implies nestedness of 
the discrete spaces, i.e. Tie C TLe+i- Moreover, we get 



H+l\T 





„l/(d-l)l, I 

q 'ne\E 



for T G T t \Tt+i, 
for T G 7/+i n Tt, 

for e g ep\ep +1 , 

for E G £p +1 n£p, 



hA E 



for e g ef\e[ +1 , 

for E G £[ +1 D£[. 



(76) 



(77) 



(78) 



We stress that 2D and 3D newest vertex bisection, e.g. [28] 
satisfy Assumption 23. 



6.2 Convergence of adaptive algorithm 

In this subsection, we show convergence of the adaptive 
scheme presented above. With the quasi-optimality (6) 
at hand, one can prove convergence of the sequence to 
a limit G "Hoc , where %oo is the closure of Ufco ^ f 
and where Uoo is the unique Galerkin solution in Hoc. 
In particular, this implies ||U^+i — U^||% — > as £ — > 00. 
A priori, it is unclear if u = U^, since adaptive mcsh- 
rcfining may lead to meshes, where the local mesh-size 
he does not tend to zero in L°°(J?) and L°°(r). 

To prove convergence, u = Uoo, we follow the esti- 
mator reduction principle of [3] as is done in [6] for (h — 
/i/2)-type error estimators and the symmetric coupling: 
By use of the Dorfler marking and the mesh-refinement 
strategy, we verify a perturbed contraction estimate 



red I 



(79) 



with certain ^-independent constants < k < 1 and 
C rc d > 0. Together with the a priori convergence, el- 
ementary calculus proves Q — Y as I — Y 00. Finally, 
reliability of the error estimator Q proves XJg — > u as 
£ — > 00. 
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One main ingredient for the proof of the estimator 
reduction estimate (79) are the following novel inverse- 
type estimates from [4] for the boundary integral opera- 
tors involved. 

Lemma 24 There exists a constant C nlv > such that 
the estimates 



\h\ /2 v rmu^n < a 



H 1 /2(r), 

h^wvtWmn <c im \\v t \\ 

||/4 /2 Vr2J^IU 2 (r) < a n v]|^]| if- 1 /2(r), 

||/lJ /2 ^ t ^|| L 2 (r) < Gnv||^|| _f/- 1 /2(_r), 

hold for all discrete functions V( £ S 1 (7e) and <P( £ 
V°(£[). The constant C- mv depends only on T and the 
7 -shape regularity of £[ . □ 

With this lemma, we can prove the main result of 
this section, which states convergence of the adaptive 
Biclak-MacCamy coupling. The same result also holds 
for the adaptive Johnson-Nedelec coupling as well as the 
adaptive symmetric coupling. The latter is treated in [4] . 

Theorem 25 Let u 6 H be the solution of the Bielak- 
MacCamy equations (30). Let XJi be the sequence of Ga- 
lerkin solutions generated by Algorithm 22. Then, the se- 
quence of corresponding error estimators Q = Pi fulfills 
the estimator reduction estimate (79). In particular, this 
implies convergence, i.e. lim ||u — U.g||-% = = lim pi. 



Proof To abbreviate notations, we define 

2- V- „.rrr\2^ V- pe{E) 2 



E 

Eern{epu£f) 



for any subset T C % U £ p U £ [ . We aim to estimate 
each contribution of the error estimator 

pI+i = Pi+i{Ti+if + P i+i(£p +1 ) 2 + pi+i(£[ +1 ) 2 - (80) 

Step 1. By use of (76), we can estimate the first 
volume contribution /9£ + i(7£+i) 2 by 



2 

Il2(t) 



\l 2 (t) 



2 

Il 2 (t) 



Pe+i(T e+ i) 2 = \\he+if\\ 2 L 2 {n} = ^ \\ht +1 f\\ 

TG71 

< E n^/iii 

TeT t nr e+ i 

+ q 2/d E UMI 

TeT e \Ti+i (81) 

= E IIWIll^T) 

-(l- 9 2 / d ) ]T iiwiii a( D 

Ter e \r e+ i 
= M%) 2 - (1 - q 2 / d ) Pe (T e \T e+1 ) 2 , 



where g denotes the mesh-reduction constant from As- 
sumption 23. 

Step 2. To estimate the second term in (80), we use 
the Young inequality (a + b) 2 < (l + 5)a 2 + (l + 6~ 1 )b 2 for 
arbitrary a, b £ M and S > 0. With this and the triangle 
inequality, we see 

pi+i{£ l+1 ) 2 



V 



<(l + 5) ^ \\K' + \[W7Ut-v]\\i HE) 
Ee£p +1 

+ (1 + <T 1 ) 



(82) 



E i 

Eesp +1 



ftjfi[(siv^ + i-av%)-i/]||^ (£) . 



Recall that 21VC/* G {P°(T e )) 2 . Therefore, [2lV£/HUi = 
[2lV£/ £ • for all sons G of a refined el- 

ement E £ £p . Furthermore, jumps over new interior 
faces vanish. With this and (77), we can estimate the 
first sum on the right-hand side of (82) by 

£ Wh^VUz ■ v]\\ 2 LHE) 

Etc (+i 



< 



E 



h] /2 pXU t -v]\\l 2(E) 



Ee£pn£p +1 
+ „!/(<*-!) 



E 



|/4 /2 [2iv[w]Hi2 (£) 



Eeep\ep +1 

= M£P) 2 -(l-q 1 /^) P nSe\£P +1 ) 2 . 

The summands in the second sum on the right-hand side 
of (82) are bounded from above by use of the pointwisc 
Lipschitz continuity of 21 and a scaling argument, i.e. 



|/»iii[(siv^ + i-avDi)-i/]||i a(B) 



(83) 



<c 2 p ||V(C/, +1 -C/,)|| 2 L2(tJf+i ^. 
We sum (83) over all interior faces and get 

Pi+i(£p+i) 2 

< (i + s)( Pe (£p) 2 (i - ? 1/(d - x) )M£iViii) 2 ) (84) 



+ {l + 5- 1 )C\\V{U e+1 -U e )\\ 



L 2 (i7), 



where the constant C > depends only on the 7-shapc 
regularity of le+i and the Lipschitz constant of 21. 

Step 3. We consider the boundary contributions of 

p 2 +1 and introduce the splitting pi{E) 2 = py(E) 2 + 

pp{E) 2 , where 



>(#) = ||^ 1/2 (0o + (i? f - - aver* ■ i/)|U>(b), 

^ a) (£) = ||^ 1/2 V r (^ -uo - ^)||l2 (£) 
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for E 6 Sf . Again we use the triangle inequality and where, as before, the hidden constant depends only on 
estimate p^\{E) by ^ an< ^ the 7-shape regularity of £[ +1 - For the last term 



in (85), we again use a scaling argument and infer 



„.i/ 2 at,* £ ll^(i(2ivc/, +1 -av^)^||! 2(s) 

+ 2\\ h t+i(*t+i-*t)h*w Z \Mu t+1 - u t )\\\ 2(n)1 

1/2, 



+ ||^ + i(2lV[/£+i - SIVC/f) ■ f|U 2 (£) where the hidden constant depends only on the Lipschitz 

fl , constant of 21 and the 7-shape regularity of Te+i- Alto- 

for E G £f + i - Summing p^^E) 2 over all boundary faces gether, we see 
and applying the Young inequality, we end up with 

< + a +7, y. iO* - <* - - »w« ■ ,)in. (E) £ (1 + "rfW - (1 - ^-■VwrvSJ") 

+ (l + S l )C(\\W(U i+1 -Ug)\\l 2(n) 
+ 3(1 + ^ 11^(^1 +||^ +1 -^||^_ 1/2(r) ) 

Bg £ r (86) 

+ 3(l + (5~ 1 ) — II /i 1 ^ 2 f^«_i_i — ^1 II 2 ' It remains to estimate the second boundary contribution 

V / / j a II £+lv £+1 t/IIL' (i?) 2a* 1 r j- 

r 4 Arguing as before, we get 

+ 3(1 + E \\hY + \(^U i+1 -^U e )-u\\l 2{E) , pfl^f 

E££f+1 (85) <(i + *) E ll^iiVr(^-«o-ar*/)||i a(B ) 

where the factor 3 stems from the inequality (E™ =1 ^) 2 < . on , x -is 111,1/2^7 (TT TT m,2 (87) 

»EJLi&i- With the help of (78), we argue similarly as + d J 2-r H ft M-i v nU/+i ~ ^JII^(B) 

before and estimate the first sum in (85) by B6£ f+i 

£ i^-ko-I)*.-"^-*,, +2(1+ ° E ? r llCVr^ (+1 -*,)|ll 1(i) . 

< V ||^ /2 (0o - - \)$t -aVU^ •^)||i 3(i 5 ) We use (78) and estimate the first sum in (87) by 

Ee£[n£f +1 

4.ai/W-D «^\^ l 2 ^ ||^ 2 Vr(^-u -^)||i 2(E) 



An inverse-type estimate from Lemma 24 can be applied Ee£[n£ r 



(88) 



to the second sum of (85). This yields 



1+1 



+ q ^p^(£[\£[ +1 ) 2 



E - ^)\\h(E } = P f\£[f (1 - g^-^pf (#V/fi) 3 - 

11/ ^ # 111 2 secon< ^ sum m (87), an inverse estimate from [7, 

= ^+1 - ^j|| i2(r) Proposition 3] can be applied, which yields 



< ||*m-^llir-v»(r) 



Here, the hidden constant depends only on r and the 
cr 



7-shapc regularity of £[ +1 - For the third sum, we use an 



E H/vdvH^+i-^lli^) 

e£[ +1 

(89) 



inverse estimate from [17, Theorem 3.6] to see = ||ftJ/_^Vr(E/£+i — Ue)\\ 2 L 2^ r ^ 

E lO^+i - ^)]] 2 o= (£ ) = lO^+i - ^)ll^(r) ~ l|c/ * +1 - Wr 



Here, the hidden constant depends only on r and the 
H ~i/2( r y 7-shape regularity of £[ +1 - Finally, the inverse estimate 
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for 23 of Lemma 24 proves 
£ ||^V r 2J(<?V+i-^)| 



2 



= \\hl£y r w($t + i-$i)\\l a w 
< \\<P e+1 -<P e \\ 2 H _ 1/2{r) 



(90) 



for the third sum in (87). Again, the hidden constant 
depends only on r and the 7-shapc regularity of £[ +1 - 
Combining (87), (88), (89), and (90) gives 

J2) (F r n2 

< (1 + <5)(pf \£[f (i - ^^)pf{e[\e^f) 

+ {l + 5- 1 )C(\\U l+1 ~U e \\ 2 m/Hr) 

+ \\<P e+1 -<P e \\ 2 H _ 1/2{r) ) 

(91) 

Step 4. We combine the estimates (81), (84), (86), 
and (91) for the different contributions of the estimator. 
This results in 

P1+1 = Pe+i(Te+i) 2 + pi + \(£ l+1 ) 2 + pe+i(£[ + i) 2 

<(i + s)(M%) 2 + Pe (£?) 2 + Pz(£f) 2 

-(l-q 1 ^)(pe(%\T e+1 ) 2 

+ Pt(£l 2 \£? +1 ) 2 + P«(£F\£f +1 ) 2 )) 
+ (l + 5- 1 )C'{\\V{U l+l -U e )\\ 2 L2{n) 
+ \\U t +i - U e \\ 2 H1/2{r) + - $4 2 H-vz(r))> 

where we have used q 2 l d < g 1 /^ -1 ), since q < 1. We 
define T t :=TiU £p U £[. Then, the norm equivalence 



/2 (n 



yields 



P 2 i+1 < (1 + - (1 - ^ d -^)p t (I i \Li +1 ) 2 ) 
+ (l + S- 1 )C\\U e+ i-Ui\\ 2 n . 



(92) 



The constant C > depends on J?, the 7-shape regular- 
ity of 7l+i and f/L l5 and the Lipschitz constant of St. 
Step 5. Recall the Dorfler marking 

< M-Mf) 2 < p e (li\l e+1 ) 2 , 

where M.t C 1^ \2^ +1 denotes the set of marked elements. 
Incorporating the last inequality in (92) gives 

^+i<(l + «)(^-(l-8 1/(,, - 1) )^) 
+ (1 + <5" 1 )C||U £+1 - U £ |||, 
= (! + *)(!- 0(1 -g^-D))^ 

+ (1 + <T 1 )C||U,+ 1 -U<|&. 



Since < q 1 /^-!) < 1 an d < 9 < 1, we have < 
l-^l-gVCd-i)) < choosing (5 > sufficiently small, 
we get < k := (1 + 5)(1 - 6(1 - q 1 ^^)) < 1, and 
Estimate (93) becomes 



< Kp 2 + C re d||U^+i — \J{ 



\h> 



with C rcd = (l + S-^C. 

Step 6. Recall that quasi-optimality (6) implies that 
linif^oo TJf G "H exists and thus ||U^+i — U^||^ —5- as 
I — > 00. Together with the estimator reduction (79), ele- 
mentary calculus predicts convergence pi — > as £ —> 00, 
cf. [3, Section 2]. The reliability ||u-U £ || H < p e of The- 
orem 12 then proves ILj — > u for £ —> 00. □ 



7 Numerical experiments 

In this section, we present some numerical 2D experi- 
ments, where we compare the three different FEM-BEM 
coupling methods on uniform and adaptively generated 
meshes. In particular, we emphasize the advantages of 
adaptive mesh-refinement compared with uniform refine- 
ment. 

Firstly, we investigate a linear problem with 21 = Id 
on an L-shapcd domain. In the second and third ex- 
periment, we choose a linear operator 21 with 21 Vm = 
(cdiff , f^) and ellipticity constant c e ii > 0. We under- 
line numerically that the assumption c c ii > 1/4 in The- 
orem 8 and Theorem 14 for the unique solvability of the 
Bielak-MacCamy and Johnson-Ncdclcc coupling is suf- 
ficient, but not necessary. Finally, we deal with a non- 
linear problem on a Z-shaped domain. Throughout, we 
consider lowest-order elements, i.e. He = Xe x with 
X { = S 1 (Tt) and y t = 7>°{8[). 

Let be a placeholder for any of the presented error 
estimators of Theorem 12, 16, or 21, i.e. Q £ {pe, 1)1 , pe}- 
The error estimator Q is split into volume and boundary 
contributions: 



ST 



a(rr=:(crr+(Qr. (94) 



res 



Recall that the variable 6 H^ 1 / 2 ^) has different 
meanings in the three coupling methods: In the Johnson- 
Nedelec and symmetric coupling </> = Vu oxt • v stands for 
the normal derivative of the exterior solution, whereas <j> 
is just a density with u ext = QJ(/> in the Bielak-MacCamy 
coupling. Quasi-optimality (6) of the coupling methods 
implies 



U/ 



H 



< 



\u - Ue\\ H i(n) 



mm 



t\\H-v 3 (r) 



(95) 



err i. 



Since the variable 4> is not comparable between the dif- 
ferent coupling strategies, we only consider the volume 
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terms and errp for comparison in the following exper- 
iments. Anyhow, we stress that one may expect that the 
finite element contribution dominates the overall conver- 
gence rate. Optimality of err^ and the corresponding es- 
timator contribution (J is numerically investigated in [5] 
for the Johnson- Nedelec coupling and some linear Lapla- 
cian in 2D. 

In the following, we plot the quantities (^P and err^ 2 
versus the number of elements N = #7f, where the se- 
quence of meshes Tt is obtained with Algorithm 22. We 
consider adaptive mesh-refinement with adaptivity pa- 
rameter 6 = 0.25 and uniform refinement, which corre- 
sponds to the case 9 = 1. For all quantities we observe 
decay rates proportional to N~ a , for some a > 0. We 
recall that a convergence rate of a = 1/2 is optimal for 
the overall error with Pl-FEM. 

Moreover, we plot the error quantities errp resp. £p 
versus the computing time tg. The time measurement is 
different for the uniform and adaptive mesh- refinement: 
In the uniform case ti, consists of the time which is used 
to refine the initial triangulation £-times, plus the time 
which is needed to build and solve the Galerkin system. 
In the adaptive case, we set i_i := 0. Then, tg consists 
of the time tg—i needed for all prior steps in the adaptive 
algorithm, plus the time needed for one adaptive step on 
the £-th mesh, i.e. steps (i)-(iv) in Algorithm 22. 

All computations were performed on a 64-Bit Linux 
work station with 32GB of RAM in Matlab (Release 
2009b). The computation of the boundary integral oper- 
ators is done with the Matlab BEM-library HILBERT, 
cf. [2]. Throughout, the discrete coupling equations are 
assembled in the Matlab sparse-format and solved with 
the Matlab backslash operator. 




Fig. 2 Volume error err £ versus numbe 
Laplace problem of Section 7.1. 
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Q BMC, unif. 






JN, adap. 
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number of elements 



N 




Fig. 1 L-shaped domain as well as initial triangulation 7o 
with #7o = 12 triangles and initial boundary mesh £q with 
#£(f = 8 boundary elements. 



Fig. 3 Volume error estimator Q versus number of elements 
N for Laplace problem of Section 7.1. 



7.1 Laplace transmission problem on L-shapcd domain 

In this experiment, we consider the linear operator 21 = 
Id. We prescribe the exact solution of (1) as 



u(x, y) = r 2 / 3 sin(|(^) 
;xt (x,y) = ilog(|x+i 



1 |2 



n 



(96) 
(97) 



on an L-shapcd domain, visualized in Figure 1. Here, 
(r, tp) denote polar coordinates. These functions are then 
used to determine the data (f,uo,4>o)- Note that Au = 
= Au cxt . We stress that u has a generic singular- 
ity at the reentrant corner. Therefore, uniform mesh- 
refinement leads to a suboptimal convergence order a = 
1/3. However, adaptive mesh-refinement recovers the op- 
timal convergence rate a = 1/2. In Figure 2 resp. 3, we 
observe optimal rates for all error and error estimator 
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3.6 - 




2.2 

10 Z 10 3 10* 10 S 10 6 10 7 

number of elements N 

Fig. 4 Efficiency index C<P/ err ? versus number of elements 
N for Laplace problem of Section 7.1. The a>axis is scaled 
logarithmically. 
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Fig. 5 Comparison of the error versus the computing time 
for adaptive and uniform refinement strategy for Laplace 
problem of Section 7.1. 



observe that the Johnson-Nedelec coupling is the fastest 
of all coupling methods, at least in this experiment. 



7.2 Linear transmission problem on L-shaped domain 

We consider a linear problem with 21 = (c n§^, ^) on 
an L-shaped domain. We again prescribe the solutions 
(w,u cxt ) by (96)-(97). Then, div2lVw = (c el] - 1)0. In 
Figure 6, we plot the error quantities eirf of the adap- 
tive schemes for different values of c e ii- We observe good 
performance of both the Bielak-MacCamy and Johnson- 
Nedelec coupling also for c c n G (0,1/4], which was ex- 
cluded by our analysis. 




10 2 10 3 10 4 10 5 10 6 

number of elements N 



Fig. 6 Convergence of the volume errors for different values 
of the ellipticity constant of the operator 21 of Section 7.2 and 
adaptive mesh-refinement. 



quantities of the different coupling methods correspond- 
ing to the adaptive scheme, whereas all quantities corre- 
sponding to uniform mesh-refinement converge with or- 
der a = 1/3. 

The quotients £jP /eirf, plotted in Figure 4, indicate 
that the residual-based estimators of Theorem 12, 16 
and 21 are not only reliable but also efficient: We see 
that (^P /err^ 2 is constant for sufficiently large N. 

When we compare the computing time, see Figure 5, 
we observe that the adaptive strategy is superior to the 
uniform one. Moreover, we see differences between the 
adaptive versions of the three coupling methods. The 
Bielak-MacCamy coupling needs significantly more com- 
puting time than the other coupling schemes. We also 



7.3 Linear transmission problem with unknown solution 

We present another experiment to underline the results 
of the aforegoing subsection. Again, we consider the op- 
erator 21 with 21 Vu = (coiiff,f|), but now on a Z- 
shaped domain, visualized in Figure 8. The data is set 
to (/, u 0} 0o ) = (1, 0, 0), and we stress that the exact so- 
lution is not known. Therefore, we plot only the error 
estimator quantities Q in Figure 9 rcsp. 10. An opti- 
mal convergence order a = 1/2 for the estimators corre- 
sponding to the adaptive schemes is observed, whereas 
uniform refinement methods lead to suboptimal conver- 
gence rates. As in Section 7.2, Figure 9 resp. 10 indi- 
cates a good performance of both the Johnson-Nedelec 
and Bielak-MacCamy coupling for c e ii G (0, 1/4]. 
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-X — BMC, c oll = 0.001 
y- BMC, c cll = 0.25 

H JN, c oll = 0.001 
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-S — sym, c cll = 0.001 
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Fig. 7 Convergence of the volume errors for different values 
of the ellipticity constant of the operator 21 of Section 7.2 and 
uniform mesh-refinement. 



Fig. 9 Convergence of the respective error estimators Q for 
different values of the ellipticity constant of the operator 21 
of Section 7.3 and adaptive mesh-refinement. 




Fig. 8 Z-shaped domain as well as initial triangulation To 
with #7o = 14 triangles and initial boundary mesh £q with 
= 10 boundary elements. 















X BMC, c oll = 0-001 
- 6 BMC, c ell = 0.25 

1 JN, c oll = 0.001 

f— JN, c cll = 0.25 
□ sym, o eU = 0.001 
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Fig. 10 Convergence of the respective error estimators Q for 
different values of the ellipticity constant of the operator 21 
of Section 7.3 and uniform mesh-refinement. 



7.4 Nonlinear experiment on Z-shaped domain 

In the last experiment we consider a Z-shaped domain, 
visualized in Figure 8, and a nonlinear operator 21 with 
2lVu = g(\\7u\)\7u, where g(t) = 2 + 1/(1 + 1) for t > 0. 
Note that the ellipticity constant of 21 is c e ii = 2. The 
prescribed solution 



u(x,y) = r 4/7 sin(|^), 

,ext/ A _ x+y+0.25 

y^iy) — (x+0.125) 2 + (y+0.125) ;i 



of (1) fulfills Z\u cxt = 0. The interior solution u has a 
generic singularity at the reentrant corner. As can be 
seen in Figure 11, resp. 12, the uniform strategy leads 



to a suboptimal convergence rate a = 2/7, whereas the 
adaptive strategy leads to the optimal convergence rate 
a = 1/2. 

The results of Figure 13 argue for the efficiency of the 
reliable error estimator Q S {pe,Tj£,Q}, which matches 
our observations of the first experiment. As in Section 7.1, 
we obtain from Figure 14 that the adaptive strategy is 
superior to the uniform one. 



(98) 

(99) 7.5 Conclusion 



In all our experiments, we observe that the adaptive 
mcsh-rcfincment strategy empirically leads to optimal 
convergence rates for the error as well as for the error 
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Fig. 11 Convergence of the volume error err e for the non- 
linear experiment of Section 7.4. 




number of elements N 



Fig. 13 Efficiency index Q /en> versus number of elements 
N for the nonlinear experiment of Section 7.4. The x-axis is 
scaled logarithmically. 
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Fig. 12 Convergence of the volume error estimator quanti- 
ties for the nonlinear experiment of Section 7.4. 
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Fig. 14 Comparison of the volume error versus the com- 
puting time for the nonlinear experiment of Section 7.4 and 
uniform resp. adaptive mesh-refinement. 



estimator quantities, whereas uniform refinement is in- 
ferior to the adaptive ones. Moreover, we see that the 
error quantities for the three different coupling meth- 
ods differ only slightly. From this point of view, one 
cannot favour a certain coupling method. But we stress 
that there are significant differences in the computing 
times, where the Johnson-Nedelec coupling is superior 
to the other two coupling methods. Although the Biclak- 
MacCamy coupling is just the transposed problem of the 
Johnson-Nedelec coupling in the linear case, it needs the 
most computing time of all three coupling schemes. Fur- 
thermore, we have numerical evidence that the reliable 
error estimators corresponding to the different coupling 
schemes are also efficient. 

Finally, we remark that our numerical experiments 
have shown that the assumption c c n > 1/4 from Theo- 



rem 8 and 14 is sufficient for the solvability of the Bielak- 
MacCamy coupling and Johnson-Nedelec coupling, but 
not necessary. 
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